Three-dimensional rogue waves in non-stationary parabolic potentials 
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Using symmetry analysis we systematically present a higher-dimensional similarity transforma- 
tion reducing the (3+l)-dimensional inhomogeneous nonlinear Schrodinger (NLS) equation with 
variable coefficients and parabolic potential to the (l+l)-dimensional NLS equation with constant 
coefficients. This transformation allows us to relate certain class of localized exact solutions of the 
(3+l)-dimensional case to the variety of solutions of integrable NLS equation of (l+l)-dimensional 
case. As an example, we illustrated our technique using two lowest order rational solutions of the 
NLS equation as seeding functions to obtain rogue wave-like solutions localized in three dimensions 
that have complicated evolution in time including interactions between two time-dependent rogue 
wave solutions. The obtained three-dimensional rogue wave-like solutions may raise the possibility 
of relative experiments and potential applications in nonlinear optics and BECs. 
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I. INTRODUCTION 



Similarity analysis is one of modern powerful tech- 
niques which allows us to find self-similar solutions 
of equations that previously were known to be non- 
integrable (see, e.g., [l| and references therein). They 
do not provide the complete intcgrability. However, 
they help to produce selected solutions in analytical 
form which may be important for variety of applica- 
tions. One of the representative examples is the nonlin- 
ear Schrodinger (NLS) equation. It is well known that in 
(l+l)-dimcnsion ((1+1)-D) this equation is completely 
integrable by inverse scattering technique 0. In (2+1)- 
D the equation is not integrable. However, solutions lo- 
calized in two transverse directions do exist Q but may 
be unstable and subjected to collapse [3, M. They are 
mostly known from numerical simulations [6[. Remark- 
ably, some of the localized solutions can be found using 
similarity reductions 0- Despite being unphysical, ex- 
act solutions provide some insight on the properties of 
the equation that is important for many applications. 
Clearly, adding a dimension changes drastically integra- 
bility properties of the equation. Thus, (3+l)-D NLS 
equation is not an exception and we are faced with the 
problem of finding its solutions knowing that they are 
not directly related to the solutions of the same equation 
in lower dimensionality. 

The NLS equation in (3+l)-D is an important model 
for variety of physical problems H, Q . It is used in non- 
linear optics Q, condensed matter physics and in par- 
ticular in modelling Bose- Einstein condensate (BEC) @ . 
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Numerical solutions can be found with various techniques 
but the value of an analytical approach is sig nificant by 
itself. In this work, extending the ideas of we use 

the similarity transformations to reduce the dimension- 
ality of the equation from (3+l)-D to (1+1)-D. In the 
former case, the coefficients in the equation are variable 
while in the latter they can be chosen to be constants. 
This allows us to use the complete intcgrability of the 
(1+1)-D equation. 

More specifically, we will focus on the possibility of 
constructing truly three-dimensional rogue waves, i.e. 
waves whose dynamics essentially depends on all spatial 
coordinates, although it is possible to identify the coordi- 
nate in which the motion is effectively one-dimensional. 
This example is directly related to the description of 
matter wave dynamics in the mean-field approximation 
(where the NLS equation is also known as the Gross- 
Pitaevskii equation), thus representing a unique possi- 
bility of creating and observing three-dimensional rogue 
matter waves. We note that the conventional rogue waves 
are either two-dimensional, as it happens, e.g. , in the 
ocean [T^], in wide aperture optical cavities jl3| and in 
capillary wave experiments [14| or one-dimensional and 
they appear in many fields including nonlinear optics (l5l - 
[L8| . cigar-shaped BECs [l9[, atmosphere [2(|, and fi- 
nances [2l| . 

The rest of this paper is organized as follows. In 
Sec. II, we describe the (l-fT)-D similarity transfor- 
mation reducing the (3+l)-D inhomogeneous nonlinear 
Schrodinger (NLS) equation with variable coefficients 
and parabolic potential to the (1+1)-D NLS equation 
with constant coefficients. In Sec. Ill, we determine the 
self-similar variables and constraints satisfied by the co- 
efficients in the (3+l)-D inhomogeneous NLS equation. 
Moreover, we give some comments about these coeffi- 
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cients. Sec. IV mainly discusses two types of localized 
3D rogue wave-like solutions, which profiles are exhib- 
ited. Finally, we give some conclusions in Sec. V. 

II. THE 3D MODEL AND SIMILARITY 
REDUCTIONS 

The original three-dimensional inhomogeneous NLS 
equation with variable coefficients can be written in a 
dimcnsionlcss form: 

= + V{r ' ^ + ^l*' 2 * + Z7 ^' (1) 

where the physical field <3> = *(r,<), r E M 3 , V = 
(dx, d y , d z ) with d x = d/dx, the external potential v(r, t) 
is a real- valued function of time and spatial coordinates, 
the nonlinear coefficient g(t) and gain/loss coefficient ^(t) 
are real-valued functions of time. This equation arises in 
many fields such as nonlinear optics (see, e.g., Q) and 
BECs (alias the three-dimensional Gross-Pitaevskii equa- 
tion with variable coefficients, see, e.g., @-Qj]]). 

We search for a similar transformation connecting so- 
lutions of Eq. (TT]) with those of the (1+1)-D standard 
NLS equation with constant coefficients, i.e. 

i^^-^^ + G^r)?^). (2) 

Here the physical field $(77, r) is a function of two vari- 
ables 77 = rj(r,t) and r = r(t) which are to be deter- 
mined, and G is a constant. Since our main goal is to 
study three-dimensional rogue waves, we choose G = — 1 
which corresponds to the attractive case (or focusing non- 
linearity in optics and negative scattering lengths in the 
BEC theory) . In order to control boundary conditions at 
infinity we impose the natural constraints (lOj 

77 —> at r-)0 and 77 —> 00 at r —> 00. (3) 

We are looking for the physical field 4'(r, t) in the form 
of the ansatz [Tl[ 

*(r, t) = ptyeW'^Qlriir, t), r(rj)] (4) 

with p(t) and <p(r,t) (like r(t), and rj(r,t), introduced 
above) being the real-value functions of the indicated 
variables, The ansatz ^ allows us to reduce the problem 
to (1+1)-D one (we notice that it differs from the one- 
dimensional stationary reductions [ijj HH). Variables in 
this reduction are to be determined from the requirement 
for the new function $(77(1-, t), r(t)) to satisfy Eq. $2$ (we 
notice that there also exist other similar reductions for 
Eq. (fTJ) which require that $(77, r) may satisfy other non- 
linear equations). Thus, we substitute the transforma- 
tion (01 into Eq. ([1) and after relatively simple algebra 
obtain the system of partial differential equations 

V 2 ri = 0, Tj t + Vip ■ V77 = 0, 2r t - |V?7| 2 = 0, (5a) 

2pt + pV 2 y - 27(i)p = 0, (5b) 

2g(t)p 2 - GIV77I 2 - 0, (5c) 

2v(r,t) + \\7^\ 2 + 2^ t =0. (5d) 



Generally speaking, equations in the system ([5]) are 
not compatible with each other when linear and nonlin- 
ear potentials are arbitrary. One, however, can pose the 
problem to find the functions v(r,t), g(t) and j(t) such 
that the system ([5]) becomes solvable. This requirement 
leads us to the procedure which can be outlined as fol- 
lows. 

• Firstly, we solve Eq. (|5"aj) subject to the boundary 
conditions ([3]) thus obtaining the similarity vari- 
ables r](r,t),T(t) and the phase <p(r,t). 

• Secondly, we consider Eqs. (|5b[) - (|5d[ ) as definitions 
for the functions p(t), v(r,t) and g(t) in terms of 
already known functions rj(r,t),T(t) and ip(r,t). 

Note that the first step determines transformation of vari- 
ables which does not involve explicitly any specific time 
dependent coefficients. However, these coefficients may 
appear after integration of these equations. The sec- 
ond step determines the coefficients which are compat- 
ible with the above change of variables. Thus, it leads 
to the model @. In this approach, the function rj(r,t) 
defines the surface where the wave has a constant am- 
plitude. The function <p(r, t) determines the wave-front 
solution (the manifold of the constant phase) . 

Thus, we can establish a correspondence between se- 
lected solutions of the (3+l)-dimcnsional inhomogeneous 
NLS equation with variable coefficients ([T]) and known so- 
lutions of completely integrable NLS equation ©. The 
latter has an infinite number of solutions thereby giv- 
ing us a chance to look for physically relevant solutions 
of the (3+l)-dimensional case. In particular, we can re- 
late them to the recently studied rogue wave solutions of 
the NLS equation [231 - 1 26l | . As a consequence, we obtain 
three-dimensional time-dependant rogue wave solutions 
of Eq. p. 

III. VARIABLES AND COEFFICIENTS OF THE 
TRANSFORMATION 

Solving Eq. (|5a[) we can write the similarity variables 
r](r,i), r{t) and the phase <p(v,t) in the form 

77(1", t) = c(t) ■ r — / c(s) ■ a.(s)ds, (6a) 
Jo 

r(t) = l£|c( S )| 2 d S , (6b) 

ip(r,t) = rf2(£)r + a(i) • r + u(t) (6c) 

where we have introduced the diagonal time-dependent 
3x3 matrix Cl(t) = dia,g(Q x (t),Q y (t),Q z (t)) with 
Cl(x(i) = — c a (t)/[2c a (t)] (hereafter a = x,y,z and 
an overdot stands for the derivative with respect to 
time). The coefficients c(t) = (c x (t),Cy(t),c z (t)), a(i) = 
(a x {t),a y (t),a z (t)) and co(t) are time-dependent func- 
tions. 
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FIG. 1: (color online). Profiles of (a) nonlinearity g(t) given by Eq. (|7b[> . (b) the coefficients A x (t) (solid line), A y (t) (dashed 
line), and A z (t) (dashed-dotted line) of second degree term of the linear potential v(r, i) given by Eq. (j8]) and (c) the gain/loss 
term 7(t) vs time for the parameters are given by Eq. (fl3|) with k x = 0.9, k y = 0.6, k z = 0.1 and k = 0.6. 



Now, from Eqs. (|5bp ~(|5d [ ) we obtain the functions 
p(£), v (r, t) and g(t) in the form 



p(t) = po ^\c x (t)c y (t)c z (t)\ e/o , (7a 
G|c(i)| 2 



2pl\c x (t)c y (t)c z {t)\e 25 ^ {s)dsl 



(7b) 



v(r, t) = v A{t) r + b(t) • r - w(t) - -|a(i)| 2 , (7c) 

where po is an integration constant. 

A few comments would be useful here. First, the 
gain/loss term j(t), is determined in the initial state- 
ment of the problem and can serve as an additional con- 
trol function or a parameter if it is a constant. Then 
changing the time-dependent dissipation we can excite 
different dynamical regimes. Second, the change of the 
all parameters is interrelated. In practical terms such a 
time dependence can be performed in different ways for 
different physical systems. In particular, in the context 
of the BEC applications, this can be done by simultane- 
ous change of the frequency of the lasers controlling the 
external trap v(r, t) and the detuning from the Feshbach 
resonance, responsible for the variation of g(t). Finally, 
we notice that one can consider the case of g(t) =const 
which however is reduced to the trivial case of a plane 
wave, whose parameters change along the chosen direc- 
tion (determined by the vector c which is a constant in 
this case). Such solutions will not be considered here. 

In writing the linear potential v(r,t) we have 
defined the diagonal time-dependent 3x3 matrix 
A(t) =diag(A x (t),A y (t),A z (t)) with the entries 



2c CT (i) cl(t) 



(8) 



and the vector function b(t) = (b x (t),b y (t),b z (t)) with 
the entries 



Ca(t) 



(9) 



It is easy to see that the velocity field v(r, t) = V<p(r, t) 
corresponding to the above-mentioned phase f(r,t) is 



given by 



v(r,t) = 2(n x (t)x, n y (t)y, ft z (t)z) + a(t) 



(10) 



such that we have the divergence of the vector field v(r, t) 
in the form 



divv(r.t) =2[n a (t) + n y (t) + n,(t)] 
= -d t ln\c x (t)c y (t)c z (t)\. 



Ill) 



Thus the zeros of any of the components of c(i) means 
the divergence of the field, which occurs at the instants 
when the nonlinearity g(t) becomes infinite [see (|7b|) ]. 
Such cases will not be considered in the present paper. 

It follows from Eqs. (|7c)) and © that if we require 
that the linear potential v(r,t) is a second degree poly- 
nomial for every space x, y, z, then we have A a ^ 0, i.e., 
c a c a — 2c^ 7^ 0, which denote that c a are not equivalent 
to constants but some functions of time. These time- 
dependent functions c a (t) will affect on the other vari- 
ables (see Eqs. (|tTa)) -r5|l) such that self-similar solutions 
of Eq. ([1]) in the form (UJ exhibit abundant structures. 
In what follows we will use specific solutions (e.g., rogue 
wave solutions) of the NLS equation to illustrate the non- 
trivial dynamics of three-dimensional rogue wave-like so- 
lutions defined by the Eq. ([T} for the different parameters 
mentioned above. 



IV. TWO TYPES OF LOCALIZED 3D ROGUE 
WAVE-LIKE SOLUTIONS 

As two representative examples, we consider the lowest 
order rational solutions of the NLS equation which serve 
as prototypes of rogue waves. First, we use the first order 
rational solution of Eq. ([2]) (see |23|). As a result, we 
obtain the first-order non-stationary rogue wave solutions 
of Eq. (JTJ) in the form 



*i(p,t) = Po ^\c x (t) c v (t) Cz (t)\eX^ ds 



1 - 



4 + 8ir(t) 



l + 2r] 2 (r,t) +4r 2 (i) 



J[ V (r,t)+r(t)] 



,(12) 
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where the variables rj(r, t), r(t) and the phase ip(r, t) are 
given by Eqs. ([5a ]) - (j6"c]) . 

For the illustrative purposes, we can choose these free 
parameters in the form 



tions of Eq. (fl) in the form 



c a (t) = a a (t) = dn(t, k a ), 

Pa = 1, 7(t) = sn(t, fc)cn(i, k), 



(13) 



(where dn, sn, and cn stand for the respective Jacobi 
elliptic functions, and k a , k are their moduli.) and w(t) = 
0. Figure [T] depicts the profiles of nonlinearity g(t) given 
by Eq. (|7b|) . the coefficients of second degree terms of the 
linear potential v (r,t) given by Eq. (|8|) and the gain/loss 
term ^(t) vs time for the chosen parameters given by 
Eq. ([TU| . The evolution of intensity distribution of the 3D 
field ([T12]) is shown in Fig. [5] We can see that the simple 
Lorenzian function of the (1+1)D case is transformed 
into a significantly more complicated evolution along the 
t-axis. The solution is localized in space and keeps the 
localization infinitely in time, which differs from the usual 
rogue wave solutions (see [? ] ) . 

On the other hand, if we choose the free parameters in 
the form 



(t)=a x {t) = l + c sin(<), 
y(t) = a y (t) = 1.2 + c cos(i), 



c. 
c 

Cz(t) 



s (t) = 0.8 + c o sin(i), c = 0.01 



(14) 



and po, Lu(t) and ^(t) are same as the ones given by 
Eq. (fT3")) , then the evolution of intensity distribution of 
the 3D rogue wave solutions (TT2]) will be changed. Fig- 
ure [3] displays the profiles of nonlinearity g(t) given by 
Eq. (|7b[) and the coefficients of second degree terms of 
the linear potential w(r, t) given by Eq. ([5J) vs time for 
the chosen parameters given by Eq. (fT4")l . For this case, 
the 3D rogue wave solution ([i"2"f is shown in Figs.|4j The 
solution is localized both in time and in space thus reveal- 
ing the usual "rogue wave" features. It is worth empha- 
sizing here, that although the "generating" function c(t) 
was chosen as monochromatic function, the respective 
change of the nonlinearity g(t) required for the existence 
of the exact solution is periodic but depending on various 
frequencies. This is natural reflection of the fact that we 
are dealing with a nonstationary solution of the nonlin- 
ear problem, characterized by the generation of multiple 
frequencies. 

Generally speaking, we have large degree of freedom in 
choosing the coefficients of transformation. As a result, 
we can describe infinitely large class of solutions of three- 
dimensional NLS equation with every exact solution of 
the one-dimensional NLS equation. Additional possibil- 
ity of choosing the solution of the latter one increases 
tremendously variety of solutions that we can obtain. 

When a hi ghe r order rational solution of the NLS equa- 
tion © (see [2J]) is applied to the transformation |QJ, we 
obtain the second-order non-stationary rogue wave solu- 



tt 2 (r,t) = P0 y/\c x (t) c y (t) c z (t)\efof(^ 

x r + P(r,,T)-iT(t)Q(r,,T) -\ e% ( r ,t)+ r( t)] (15) 
H(t),t) J 

where these functions P(?7,t), Q(?7,t) and H(tj,t) are 
given by [24[ 



P(t],t) = -\rf - 6t] 2 t 2 - 10t 4 - ^if - 9r 2 



Q(V, r) = rf + ir) t + 4t 4 - 3rf + 2t a 
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(16) 



H{ V , r) = ^rf + \fr 2 + r, V + + \rf 

9 a 3 o o 9 n 33 o 3 

H — T n T H 77 H T H . 

2 2 ' 16 ' 8 32 

The variables rj(r,t), r(t) and the phase tp(r,t) here are 
given by Eqs. ([paj) -(|5c l) . 

As in the previous cases, we choose the parameters 
given by Eqs. (TJ3]) and ((HI) except for a a (t) = 0. The 
intensity distributions of the second-order rogue wave so- 
lutions (TT51) arc depicted in Figs. [S] and [51 Clearly, the 
field evolution in this case is more complicated. In one 
case, shown in Fig. [51 the solution is localized in all three 
dimensions in space. In the other case, shown in Fig. [6] 
the solution is localized both in space and in time thus 
displaying the basic feature of a rogue wave that "appears 
from nowhere and disappears without a trace" . 

It follows from the above-mentioned two cases for the 
parameters that the parameters c a (t), a CT (t) and ^(t) 
can be used to control the wave propagations related 
to rogue waves, which may raise the possibility of rela- 
tive experiments and potential applications in nonlinear 
optics and BECs. Similarly we can also obtain three- 
dimensional higher-order time-dependent rogue wave so- 
lutions of Eq. ((TJ) in terms of the transformation ((4} and 
higher-order rogue wave solutions of the NLS equation 
@, which are omitted here. 

As always happen with the nonlinear Schodinger equa- 
tion in two and three dimensions, their localized solutions 
may collapse. Stability of the solutions presented in our 
work is still an open question. This question deserves 
separate studies as it is a task that is far from being 
trivial. We leave these studies to later publications. 



V. CONCLUSIONS 

In conclusion, we have presented similarity reduc- 
tions of the (3+l)-dimcnsional inhomogeneous nonlinear 
Schrodinger equation with variable coefficients to (1+1)- 
dimensional one with constant coefficients. This trans- 
formation allows us to relate certain class of localized so- 
lutions of the (3+l)-dimensional case to the variety of so- 
lutions of integrablc NLS equation of (l+l)-dimcnsional 
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FIG. 2: (color online). Color coded plot of wave intensity (a) |*i \ 2 (x, 0, 0, i) with max {Xj0?0>t } |*i| 2 = 0.06, (b) *i j 2 (0, y, 0, t) 
with max{ 0i!/i o,t} |*i| 2 = 0.068 and (c) |*i| 2 (0, 0, z, t) with max {0i0iZ . t } ]*i| 2 = 0.068 defined by the solution (JT5J for the 
parameters given by Eq. ()13[l with k x = 0.9, k y — 0.6, k z = 0.1, k = 0.6. Note that for given set of parameters, the wave is 
localized in space. 
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FIG. 3: (color online). Profiles of (a) nonlinearity g(t) given 
by Eq. l[fb]) and (b) the coefficients A x {t) (solid line), A y (t) 
(dashed line), and A z (t) (dashed-dotted line) of second degree 
term of the linear potential v(r, t) given by Eq. (j8]) vs time 
for the parameters are given by Eq. (|14|l . 



case. As an example, we illustrated our technique by 



two lowest order rational solutions of the NLSE. These 
are transformed into rogue wave solutions localized in 3D 
space that have complicated evolution in time. The tech- 
nique may also be extended to other NLS-type equations 
to exhibit their rogue wave solutions. 
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